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We propose an approach for single spin measurement. Our method uses techniques from the 
theory of quantum cellular automata to correlate a large amount of ancillary spins to the one to 
be measured. It has the distinct advantage of being efficient, and to a certain extent fault-tolerant. 
Under ideal conditions, it requires the application of only O(viV) steps (each requiring a constant 
number of rf pulses) to create a system of N correlated spins. It is also fairly robust against pulse 
errors, imperfect initial polarization of the ancilla spin system, and does not rely on entanglement. 
We study the scalability of our scheme through numerical simulation. 
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One of the most interesting challenges in physics to- 
day is that of measuring the state of a single (nuclear) 
spin. Being able to do this would bring us closer to spin 
based quantum computers , and have a myriad of ap- 
plications, ranging from spintronics, to protein analysis. 
Unfortunately, performing such measurement is not an 
easy task. Several methods have been proposed 0, and 
single-spin detection has been done |3( , and the measure- 
ment in some specific cases has also been achieved 

Cappellaro et. al. 5] propose using a system of N 
independent spins as a measurement system. This sys- 
tem is coupled to the spin being measured, and using en- 
tangling operations, creates a large correlated state, one 
whose signal can be measured with current NMR tech- 
nology. The methods presented there require O(N) pulse 
sequences in order to achieve N quanta of polarization. 
We expand on their ideas to create a scheme which uses 
0{^/~N) pulses, and is fairly robust against noise and er- 
rors, and does not rely on entanglement. The improved 
running time is the most important advantage, since the 
whole procedure must finish before decoherence destroys 
the information being measured. 

The method presented here is inspired by quantum cel- 
lular automata 0, 0] and pulse driven quantum comput- 
ers H, 0- It uses a cubic lattice crystal 01 with two 
nucleus types, which we call A and B. We assume, for 
the time being, only nearest neighbor interaction. Each 
species A nucleus is connected only to B nuclei, and vice 
versa in a checkerboard fashion. See fig. ^ 

We will refer to upward and downward z polarizations 
as | + 1) and | — 1). We assume, for the time being, that 
the crystal is initialized to a completely polarized state 
with all nuclear spins in a downwards z polarization state, 
| — 1). The method then consists of bringing one corner 
of the crystal into close proximity to the spin we wish to 
measure, so as to couple the two spins. Once coupled we 
can use NMR rf pulses to correlate them, or swap the 
states. Suppose the spin we wish to measure is initially 
in the state \ip) = a\ — 1) + j3\ + 1). After the swap the 
top-left vertex nuclear spin will be in the state \ip). 

What we present now is an efficient method to create 
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FIG. 1: Cube lattice: A crystal with two types of nuclei A 
and B, one represented as light gray spheres, the other dark 
gray. Each species A is neighbored by only B type nuclei 
and vice-versa. The lines connecting the spheres represent 
the nearest neighbor couplings. The white sphere represents 
the spin we wish to measure. This spin is coupled to the dark 
gray nucleus in the top-left vertex. 



a very large correlated state within the crystal. Under 
ideal conditions, that is complete polarization, perfect 
pulses, and no decoherence, the method creates the state 
a\ - 1)® N + j3\ + 1)® N using only 0{yfN) steps (each 
requiring a constant number of rf pulses). When N ~ 10 6 
the resulting state gives a strong enough magnetic signal 
to be measured. Achieving this polarization, in the ideal 
case, requires applying about 200 steps of our algorithm. 

In order to more easily illustrate our algorithm it is 
best to visualize the cube lattice in the following way. We 
envision slicing the cube into layers, such that the first 
layer is the corner nuclear spin that contains the state 
to be measured. Layer two contains all nuclei coupled 
to layer one. Layer three contains all nuclei coupled to 
layer two which are not in layer one, and so on. This is 
illustrated in fig. |21 

Each layer includes spins of only one species, layer one 
being all A, layer two all B, layer three all A again and 
so on. We envision taking only half the cube lattice, so 
that each layer is larger than the previous one. Layer i 
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FIG. 2: Pyramid lattice: A different view on the same struc- 
ture as fig. Q The top gray sphere is layer 1 and corresponds 
to the top-left vertex spin in fig. The three light gray 
spheres directly below it are layer two, and so on. The lines 
connecting the spheres represent the nearest neighbor cou- 
plings, and are exactly the same as in fig. 



has i more spins than layer i — 1. 

For sake of analysis, suppose the crystal's Hamiltonian 
has nearest neighbor (Manhattan distance one)couplings 
only. In such a Hamiltonian, each A spin is coupled only 
to B spins, and vice versa. The resonant frequency of 
each A (conversely B) spin is affected by the states it's 
B (conversely A) spin neighbors. Since the neighbors 
are all indistinguishable spins, only the total field value 
is important, i.e. the number of neighbors pointing up, 
minus the number of neighbors pointing down. 

Using the methods first suggested by Lloyd also 
used in j§| and others, it is possible to create gates that 
address species A spins that have a particular neighbor 
field value. For instance, it is possible to apply a NOT 
gate (using the language of quantum computation) to all 
A spins that have neighbor field value 0, and only these 
spins. 

We now show how to use the lattice structure to 'am- 
plify ' the spin we wish to measure to a detectable signal 
strength. We call the set of all species X spins, with 
neighbor field fc, X/,. 

The algorithm is the following. Repeatedly apply a 
NOT gate to the following sets of spins: -B-2, B-x, Bq, 
A-2, A-2, Aq. Each NOT gate consists of a (constant 
size) sequence of strongly modulated pulses. 

Suppose that a = 1, i.e. \tp) = \ — 1). Then absolutely 
nothing happens to the lattice (ignoring errors), and so 
it remains in the all | — 1) state. To see this, note that 
all lattice points have at least three neighbors, all in the 
| — 1) state. Hence the field k of every cell is at most —3. 
Since we are only doing flips on k = 0,-1,-2 this does 
not affect the lattice at all. 

Suppose now that = 1, i.e. the top cell of the pyra- 
mid is initialized to the state | + 1). Then, all B neighbors 
of this vertex (those in layer 2) will have field value —2: 
—3 from three downward neighbors each in state | — 1) 
and a +1 from the upward neighbor in state | + 1) . There- 
fore they will be flipped on the first it pulse. No other 
nuclear spins will be affected. In the next A stage, all A 
cells in layer 3 will be flipped, and so on. In n stages all 
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FIG. 3: Ideal NMR Spectrum: This is the ideal spectrum for 
a spin on the second layer. There are 5 distinct peaks, one 
for each of the five neighbor field values: —4, —2, 0, +2, and 
+4. 



cells in the first n layers will be flipped to | + 1). 

After n stages there will be N = h(n + l)n(n — 1) 

nuclear spins pointing up. Hence, only 0(\^N) stages 
are required to obtain a total field of N spins. 

By linearity, in the general case where \ip) = a\ — 1) + 
f3\ + 1), in the complete absence of decoherence the pro- 
cedure would result in a cat state of N nuclear spins in 
the state a\ - 1) 8JV + /3\ + 1)® N . 

At the end of such procedure, one can measure the 
magnetic field of the cube lattice crystal, obtaining the 
desired measurement result. 

Khitrin et. al. used a similar approach to ours in 
their scheme for polarizing spin chains. Although their 
method could potentially be used for spin measurement, 
the method presented here is cubically more efficient. 

In order for the scheme to work it is tantamount that 
the frequencies of the target spins be different from the 
ones we do not wish to affect. It would also be highly 
desirable to have the resonant frequencies of all target 
spins lumped together, as opposed to interspersed with 
the non-target frequencies. In the ideal setting, where 
the Hamiltonian has only first neighbor (strictly speaking 
Manhattan distance 1) couplings, this is the case. In 
figure|31we see what an ideal spectrum for a second layer 
spin looks like. It has 5 distinct peaks, one for each of 
it's possible neighbor fields: —4, —2, 0, +2, and +4. 

We have so far assumed a nearest neighbor coupling 
Hamiltonian. A solid crystal, however, is governed by a 
direct dipole coupling Hamiltonian, 



where fcjj equals one if i, j are of the same species, and 
zero otherwise. And 
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where r^j is the distance between the two nuclei, @ij is 
the angle between the vector connecting the two nuclei 
and the z-axis (determined by the magnetic field) and 
gij is a simple constant that depends only on the nuclear 
types of i and j . 
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FIG. 4: Rhombohedral lattice NMR Spectrum: This is the 
spectrum for a spin on the second layer of a rhombohedral 
lattice with angles a = j3 = 7 = n/3. The non-nearest neigh- 
bors, and in particular, the homonuclear couplings, swamp 
out the peaks. It is impossible to discern and address the 
wanted spins. 




FIG. 5: Suppressing homonuclear coupling: This is the spec- 
trum for a spin on the second layer of a rhombohedral lattice 
with angles a = = 7 = 7r/3, with the homonuclear cou- 
plings suppressed. Note the similarity to the ideal case. 

Since the angle between the edges of the lattice and 
the z-field affect the coupling strength, it makes sense to 
orient the crystal such that the top of the pyramid points 
up in the z direction, making the angles symmetric. The 
problem is that, in a cubic lattice (with the Bravais angles 
a = (3 = 7 = tt/2), all edges in the lattice have an angle 
to the z pole of exactly arctan \/2. This is known as the 
magic angle. At this angle the coupling becomes zero. 

It is clear that a cubic lattice is inappropriate if we 
are in a dipole coupling regime, and we wish to have the 
angles symmetric. A simple solution is to use a rhombo- 
hedral crystal structure. We take one where the Bravais 
angles are a — ft = 7 = 7r/3. In figure 0] we show the 
spectrum of a second layer spin. Same as before, the 
pyramid peak points up in the z direction. Although the 
nearest coupling are now there, they are washed away 
by the homonuclear couplings and non-nearest neighbor 
couplings. It is clear that we cannot address the spins we 
wish. 

However, it is possible to suppress all homonuclear cou- 
plings(c.f. 11]). When we do so we get a Hamiltonian 
that looks remarkably close to the ideal nearest neighbor 
one. In figure [5] we show the spectrum of a second layer 
spin. Notice how similar it is to the ideal case. We can 
show similar results for different Bravais angles; in con- 
clusion, almost any crystal structure will work properly 
when we suppress homonuclear interactions. 

In all cases presented, the crystals are hypothetical 
constructs, where the couplings are set up so that the 



nearest neighbor couplings are roughly 1000 Hz. 

It is clear that imperfect polarization will reduce the 
efficiency of our algorithm, yet it is a necessary considera- 
tion. Also, we have so far assumed only nearest neighbor 
interaction, which is not usually the case. A weak cou- 
pling can be observed between spins at greater distances, 
not taking these into account will result in imperfections 
in our gates. It is important to deal with these issues. 

One important fact to notice is that our scheme is im- 
pervious to phase-flip errors in the lattice. Effectively, all 
pulses, and the final output, depend only on the diago- 
nal terms of the density operator p of the lattice. This 
is critical to state, not only because it deals with the al- 
gorithm's robustness to errors, but also because it allows 
us to do an approximate, non-coherent, simulation of the 
evolution of our algorithm, as we show below. 

The diagonal elements of p must be protected, in the 
language of quantum information this amounts to pro- 
tecting p against bit-flip errors. A single bit-flip error in 
the cube lattice by itself does not drastically change the 
overall field, however, if this bit flips during the running 
of the algorithm it has the potential to influence all lat- 
tice points in a radius around it during the next step of 
the algorithm, by changing their neighbor fields during 
the pulse sequence, potentially creating a cascading of 
errors effect. 

The first fact to notice is that errors inside the pyramid 
are less critical than errors on the faces, or worse, edges, 
of the pyramid. The reason is that inside the pyramid 
all lattice points have exactly 6 neighbors. If one flips 
erroneously its neighbors still have 5 neighbors each in 
the correct state. 

Take an interior lattice point that is in the state | — 1) 
and should be flipped to | + 1) in the next pulse sequence. 
Since it is supposed to flip it should have three positive 
neighbors and three negative ones. Suppose that an up- 
ward neighbor has an error, then it will have two positive 
and four negative. This is, however, not a problem. The 
lattice point in question will have a field value k = — 2, 
instead of k — 0, but will still be flipped correctly. 

Contrast this to the case where a bit flip occurs along 
the edge of the pyramid. In this case the field of the 
lower neighbors to a lattice point with error will go from 
having a field of —2 to a field of —4 and will no longer 
be correctly flipped at all. The error will then propagate 
downwards. 

One way to reduce errors is to extend the flipping op- 
eration to spins with neighbor fields of +1. This will 
correct some edge-related errors. One other way to dras- 
tically reduce error propagation is to force errors into the 
inside of the pyramid. 

Previously we concluded that we need to suppress the 
homonuclear coupling during our pulse sequences. In this 
Hamiltonian, the term a + ' a_ + a_ VV , often referred 
to as the 'flip-flop ' term, has the effect of swapping two 
anti- aligned neighboring spins. Take a horizontal slice of 
the pyramid and suppose there is a single error among 
these lattice points. The effect of the flip-flop term on 
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FIG. 6: Output Signal: This graph shows the signal strength 
at the end of the algorithm, as a function of the lattice size. 
The size is given in number of nuclei, and the signal strength 
is in number of correctly polarized nuclei. These numbers are 
for an initial crystal polarization of 90% and gate error rates 
one and five percent. 

neighboring spins is well known and is called spin diffu- 
sion, |12|. A first order approximation to this dynamic 
is a quantum walk on the lattice. In this model, an er- 
roneous bit not only becomes diffused over several spins. 
Also, it has a much lower probability (averaged over time) 
of being measured in non-interior points. This leads us to 
believe that letting the system evolve under the homonu- 
clear coupling for some time before running the algorithm 
will have the effect of reducing errors in the edges, and 
increase overall robustness to errors. 

In order to test the reliability of our algorithm in the 
large scale it would be infeasible to have a coherent sim- 
ulation. However, as stated above, the non-coherent na- 
ture of our procedure allows us to take a different ap- 
proach. One simulation we set up is a non-coherent 



Monte Carlo style simulation. We set up a three dimen- 
sional byte array, each byte representing a nuclear spin. 
Each byte is set to an initial value of —1. In order to 
simulate an initial polarization of 1 — eo we then flip the 
sign of each byte with probability eo. Each pulse is simu- 
lated in a similar fashion. A gate error ei is simulated by 
failing to flip bytes with probability e\. This simulation 
is repeated a number of times in order to get a meaning- 
ful statistical average. This technique allows to simulate 
very large systems, of size up to 10 s nuclear spins. 

We calculated the signal strength at the end of our 
algorithm under different conditions of initial crystal po- 
larization, and pulse error rate. In figure we plot the 
signal strength as a function of lattice size, when the ini- 
tial polarization is 90% and pulse error rates are 1% and 
5%. As can be seen from the graph, in order to get a 
signal strength in the order of 10 6 polarized spins it is 
necessary to use a crystal lattice with about 10 8 nuclei. 
Running our algorithm on such a lattice requires about 
800 gates. 

In closing, moving to a three dimensional measurement 
system, as opposed to a ID chain as in [5j, gives distinct 
advantages. The clearest benefit is a cubic speedup in 
the procedure. This speedup was expected. Also the 
algorithm gains substantial robustness to errors. This 
is an even greater benefit, given the current fidelity of 
NMR operations. Our numerical simulations tell us that 
the scheme is highly efficient and robust. However, the 
nature of our simulations do not allow us to take into 
account large coherent superpositions or errors through- 
out the pyramid lattice. Large scale coherent simulations 
are infeasible. Ultimately we will need an actual physical 
implementation in order to best understand the capabil- 
ities, and limits, of our proposed scheme. 

This paper was written, in part, with funding from 
ARDA, ARO, CFI, CIAR, DARPA, LPS, MITACS, 
NSERC, and ORDCF. 
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